Forward-backward elliptic anisotropy correlation in parton cascade 
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A potential experimental probe, forward-backward elliptic anisotropy correlation (Cfb), has been 
proposed by Liao and Koch to distinguish the jet and true elliptic flow contribution to the measured 
elliptic flow (V2) in relativistic heavy-ion collisions. Jet and flow fluctuation contribution to elliptic 
flow is investigated within the framework of a multi-phase transport model using the Cfb probe. 
We found that the Cfb correlation is remarkably different and is about two times of that proposed 
by Liao and Koch. It originates from the correlation between fluctuation of forward and backward 
elliptic flow at low transverse momentum, which is mainly due to the initial correlation between 
fluctuation of forward and backward eccentricity. This results in an amendment of the Cfb by 
a term related to the correlation between fluctuation of forward and backward elliptic flow. Our 
results suggest that a suitable rapidity gap for Cfb correlation studies should be around ± 3.5. 
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The results from Brookhaven Relativistic Heavy-Ion 
Collider (RHIC) indicate that a strongly-interacting par- 
tonic matter has been created in relativistic nucleus- 
nucleus collisions [l| . Two powerful probes exposing the 
characteristics of the new matter are elliptic flow and jet. 
Elliptic flow, has been measured via the second Fourier 
coefficient (i>2) in the azimuthal distribution of final par- 
ticles [2, 3] . The V2 data show remarkable hydrodynam- 
ical behaviors, which implies the formed matter is ther- 
malized in a very short time and expands collectively as 
a liquid with low shear viscosity/entropy. On the other 
hand, jet, which is produced non-collectively by initial 
hard scatterings, has been experimentally studied by nu- 
clear modification factor and jet-like correlation |4I-[l0|. 
These observation shows that jet losses energy when it 
passes through the hot and dense QCD medium formed 
in heavy-ion collisions. However, two probes, elliptic flow 
and jet, are correlated in fact. The jet contributes to the 
anisotropy in final azimuthal distribution and certainly 
to the observed v-i , especially for high transverse momen- 
tum (pr) because of the path length dependence of en- 
ergy loss in non-central collisions [ll|, [l2| ■ Nuclear com- 
munity usually classify such a non-collective behavior as 
non-flow. Therefore, it is important for elliptic flow mea- 
surement to distinguish these two parts from observed V2 
to further separate collective and non-collective proper- 
ties of the new matter. A potential experimental probe, 
forward-backward elliptic anisotropy correlation (Cfb), 
has been proposed for this purpose by Liao and Koch 
(LK) [l3| . It is believed that the new probe can distin- 
guish how jet and elliptic flow contribute to the observed 
V2 which n probably reveal the momentum scale where 
elliptic flow or (semi-) hard processes dominate. Using 
a two-component parameterization, they found the F-B 
correlation (Cfb), which takes a maximum value of unity 



at low transverse momentum (pr) and falls to zero with 
Pt- This indicates the jet contribution to the measured 
i>2 grows with increasing px and even jet contribution 
dominates at high px- However, event- by- event fluctua- 
tion of V2 is also important for v 2 measurement [l2l [l4| . 
Experimental data shows that the relative nonstatistical 
fluctuations of the V2 is to be approximately 40% 14 1. 
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Unfortunately it is not taken into account in LK's cal- 
culation. Our present study considers the effect of V2 
fluctuation to further understand bulk properties and 
related jet effects in the early stage of heavy-ion colli- 
sions. Our paper presents the relative nonstatistical fluc- 
tuation effects on the forward-backward V2 correlation 
(Cfb) and the G- factor range which reflects jet contri- 
bution to observed V2 with a multi-phase transport model 
(AMPT) [lj| . A more generalized Cfb which includes 
F-B elliptic flow fluctuation correlation is formulated. 

The AMPT model consists of four main components: 
the initial conditions, partonic interactions, conversion 
from partonic to hadronic matter and hadronic interac- 
tions. The initial conditions, which include the spatial 
and momentum distributions of minijet partons and soft 
string excitations, are obtained from the Heavy Ion Jet 
Interaction Generator model. Scatterings among partons 
arc modeled by Zhang's Parton Cascade model, which 
at present includes only two-body scatterings in which 
cross sections are obtained from the pQCD calculation 
by screening mass. In the version with string melting 
mechanism, partons include minijet partons and partons 
from melted strings. The quark coalescence model is used 
to converte partons into hadrons. The dynamics of the 
subsequent hadronic matter is then described by a rela- 
tivistic transport model. The detail of the AMPT model 
can be found in Ref . 15[ . From the previous AMPT cal- 
culations, it is found that elliptic flow can be built by 
strong parton cascade @, [l5|-[lj| and jet losses energy 
into partonic medium to excite a Mach-like cone struc- 
ture (2(jl - f22l ]. It is clear from the above that partonic 
effect can not be neglected. Therefore the string melting 
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AMPT version is appropriate when the energy density is 
much higher than the critical density which is predicted 
for phase transition. In this work, we use the string melt- 
ing AMPT model to simulate Au+Au collisions at ^/sjvjv 
= 200 GeV. The partonic interaction cross section is set 
to 10 mb. 

Elliptic flow v 2 is defined as 



We introduce a G factor to represent total jet contribu- 
tion V 2 to total V2(pt), which has a relation to g factor 
in llJl as followed: 



G = 



> (7) 

v 2 



Jo 27r #cos20 < > _< V 2 (p t ) > 



< v 2 (pt) > 



IT d< t> < 



d 2 N 
ptdptdtp 



< my > 
> *- p t dpt ^ 

(1) 

where the total yields d , N ,, arc resulted from both cl- 
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liptic flow and jet |l3[ which then contribute to the mea- 
sured v 2 . 

The proposed observable Cfb[pt] is the correlation of 
the total elliptic flow V 2 [prjbetween forward (F) and 
backward (B) rapidity bins pjj|, defined as: 



Using a different form of total V 2 (j>t) in forward (F) 
and backward (B) rapidity bins, i.e. V F = (| +r/ F )(l — 
G)(V 2 ) +ZG(V 2 ) and V B = (| + V B )(1 - G)(V 2 ) + (1 - 
£,)G(V 2 ), we get the formula below with the assumption 
(if) = ( V B ) = 0, (0 = (1 - = |, and C(l -0=0, 
but (?7 F ?7 B ) ^ 0: 
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With the assumption (rj F ) = (r/ B ) = 0, (£) 
= 5 an d £(1 — — 0' changes to 
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where r} F ^ represents random deviations from the aver- 
age elliptic flow yield in forward (backward) rapidity bin, 
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When the correlation of the fluctuation of F-B elliptic 
flow is zero, (i.e. (rj F ri B ) = 0), C FB = 1 - G 2 . From 
here, only jet will contribute to v 2 when G = 1, Cfb = 0. 
On the other hand, only elliptic flow contributes to v 2 
when G = 0, Cfb = F While when the correlation 
is nonzero, for example (?7 F ?/ S ) = \, Cfb = 2(1 — G). 
Then we can find that only jet contributes to v 2 with a 
condition G = 1, Cfb = 0. However, only elliptic flow 
will contribute to v 2 when G = 0, Cfb = 2. 

The current «2 measurement can not distinguish the 
deviations of elliptic flow rj F ^ from jet effect of £ in ex- 
periment, since they are mixed together. One can rede- 
fine total V 2 (pt) in forward (F) and backward (B) rapid- 
ity bins, V F = (l + vf)(V 2 F ) , and V B = (l + <)(y 2 s ), 
where "HfIb) includes the fluctuation of elliptic flow and 



represents elliptic flow (jet) part in jj 2 , ( represents J et contribution. The F-B correlation can be rewritten 



the jet contribution to F or B rapidity bin, and g factor, 
giving the relative weight of the jet contribution to the 
total, i.e. (F+B), yield, is defined as [lU 
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When the correlation of the fluctuation of F-B elliptic 
flow yield is zero, (i.e. (rj F rj B ) = 0), then LK found: 
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When only jet contributes to v 2 , that is g — 1, Cfb = 0; 
On the other hand, if only elliptic flow contributes to v 2 , 
then g = 0, Cfb = F 

However, if the correlation of the fluctuation of F-B 
elliptic flow is nonzero and for example (i] F i] B ) = \, 
Cfb will change to 
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When only jet contributes to v 2 , then g = 1, Cfb = 0; 
However, if only elliptic flow contributes to v 2 , then g = 
0, Cfb = 2. 



With the assumption (r] F 11 ) = (r) B u ) = 0, we can finally 
get the same expression as Eq. ©. 

Cfb has been simulated by using AMPT model for 
Au+Au at v'snn = 200 GeV at the centrality bin 0-20%. 
The results are shown in Fig. [I] where different symbols 
represent different forward and backward pseudo-rapidity 
bins. For example, 0.5 < |t7| < 1.0 represents that the 
forward bin is within 0.5 < r\ < 1.0 and while the back- 
ward bin is within —1.0 < r\ < —0.5. Cfb remains 
constant for different rapidity gaps at low px (j>t < 2 
GeV/c) in 0-20% centrality bin. To our surprise, Cfb 
appears to be remarkably different from the prediction 
of Ref. [l3|], nearly two times as LK's (unity) at low pt 
(pt < 2 GeV/c) for the mid-rapidity gap, and decreases 
with rapidity gap. The reason, as shown in Eq. ([5]) and 
©, is that the correlation of the fluctuation of F-B el- 
liptic flow (rj F ri B ) is nonzero, i.e. (rj Fr q B ) is close to | 
at low px for mid-rapidity gap at 0-20% centrality bin. 
As rapidity gap rises to 7.0 (Fig. Eh)), Cfb falls to 1.2, 
which indicates (j] F J] B ) becomes weak for a large rapidity 
gap. 
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FIG. 1: Cfb for 200 GeV/c Au + Au collisions at 0-20% rapidity gap 0.5 < |r/| < 1.0 at centrality bin 0-20%. 
(circles), 20-40% (squares), and 40-60% (triangle) centrality 
bins, (a) 0.5 < \q\ < 1.0, and (b) 3.0< \q\ <4.0 (Some symbols 
are slightly shifted in pT-axis for clarity). 



Fig. Q] also displays that Cfb is much larger for cen- 
tral collisions than non-central collisions at low px, which 
indicates that F-B elliptic flow fluctuation is strongly cor- 
related in the most central collisions. This trend is con- 
sistent with centrality dependence of elliptic flow fluc- 
tuation 24 1 . In addition, (f] F r] B ) is near to a positive 



constant at low px for the given rapidity bins and cen- 
trality. However, at high px {px > 4 GcV/c), Cfb is 
close to zero, indicating that only jet contributes to the 
observable «2 and (r] F rj B ) is zero. The results are similar 
with LK's finding here. Therefore, (r] F r) B ) should posi- 
tively decreases in intermediate px region (2 < px < 4 
GeV/c). It means that Cfb can not be used to directly 
measure the jet contribution to observed V2 without the 
knowledge of (r] F r] B ), especially for mid-rapidity gap. 

The same tend is also seen in Fig. [Tfb) , rapidity bin 
3.0 < \r)\ < 4.0. For centrality bins 20-40% and 40-60%, 
Cfb looks near to unity at low px, which means (rj F rj B ) 
is close to zero. The structures are similar with LK's 
finding here. Therefore, Cfb with large rapidity gap is 
a clean probe to extract jet contribution to elliptic flow, 
as Ref.[l3( expected. 

As we discussed in previous section, two limits of Cfb, 
i.e. Cfb = 1 — G 2 and Cfb = 2(1 — G) which can be 
deduced assuming (rj F r] B ) approaching to zero at high 
px or (i] F ri B ) approaching to j at low p T , respectively. 
Fig. [2] shows the above two limits. We can see that the 
factor G, which reflects the jet contribution to the ob- 
served V2, cannot be strictly quantified by the observed 
Cfb- This gives only a range of G- value which can be 
estimated for mid-rapidity gap. This is due to the reason 
that we do not know the quantitative relation of (r) F r) B ) 
and G for intermediate px region. With the decreasing 
of Cfb, the uncertain range of G becomes narrower. 

On the other hand, elliptic flow, as is shown in 
Ref . 0, H|| , is converted from initial collision geometry. 
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FIG. 3: Cfb for initial eccentricity (open symbols) and el- 
liptic flow (solid symbols) in Au + Au collisions at y / ijvjv = 
200 GeV at centrality bin 0-20% within different rapidity gaps 
represented by different symbols. Some symbols are slightly 
shifted in pr-axis for clarity). 



It means that fluctuation of elliptic flow should reflect 
the information on fluctuation in the initial state geome- 
try. Moreover, the fluctuation due to the magnitude and 
centrality dependence of observed elliptic flow are con- 
sistent with fluctuation of the initial geometry shape of 
the collision region (26|. It has been found that such a 
initial geometry irregularity can be transferred into fi- 
nal momentum anisotropics by strong partonic interac- 
tions [H, (eg. triangular flow). It becomes more in- 
teresting to compare Cfb for elliptic flow (112) and initial 
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eccentricity (e2). 

In AMPT model, partons, which are melted from 
string, are initially distributed in momentum and coor- 
dinate spaces. We can similarly define Cfb for initial 

eccentricity (e 2 ) as Cfb^) = ' wnere E 2 ^ 

are the total eccentricity (£2) in forward (F) and back- 
ward (B) rapidity bins, and (E 2 (p T )) = <g^(-™-}. 

Similar to V 2 (B) , E 2 [B) can be written as E 2 = 
(5 + Vf)( e 2) and Eg = (| + ri B ){E 2 ), where r^ (B) 
represents random deviations from the total eccentric- 
ity in forward (backward) rapidity bin. If one compare 
with the former analysis of v 2 , Cfb{c2) presents similar 
trend to Cfb (1*2) , if the correlation between fluctuation 
of forward and backward elliptic flow is mainly due to 
the initial correlation between fluctuation of forward and 
backward eccentricity. 

As shown in Fig. [31(a), it is an obvious fact that 
Cfb{v 2 ) is efficiently transferred from Cfb^) at small 
rapidity gap. However, at large rapidity gap (in 
Fig. [3(b)), Cfb(v2) is suppressed much deeply in con- 
trast to CFB(e 2 ). This indicates that parton cascade are 
not strong enough for larger rapidity gap to transfer the 
F-B correlation of initial e 2 fluctuation to final F-B cor- 
relation of v 2 fluctuation. This is because of the fact that 
parton interactions are weak at higher rapidity. There- 



fore, the measurement of Cfb for elliptic flow (v 2 ) may 
give us more information about the correlation of initial 
geometry fluctuation. 

In conclusion, jet and flow fluctuation contribution to 
elliptic flow are investigated by forward-backward ellip- 
tic anisotropy correlation (Cfb) within the framework of 
a multi-phase transport model. We found that the Cfb 
correlation is remarkably different from LK's and is about 
nearly double of that proposed by Liao and Koch. It 
stems from the correlation between fluctuation of forward 
and backward elliptic flow at low transverse momentum, 
which is mainly due to the initial correlation between 
fluctuation of forward and backward eccentricity. This 
leads to an amendment of the Cfb by a term related to 
the correlation between the fluctuation of forward and 
backward elliptic flow. The present study shows that 
Cfb correlation decreases with rapidity gap and a suit- 
able rapidity gap for Cfb correlation studies should be 
around ±3.5. 
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